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Abstract 

This paper is concerned with the study of one-body dissipation effects in 
idealized models resembling a nucleus. In particular, we study the quantum 
mechanics of a free particle that collides elastically with the slowly moving 
walls of a Bunimovich stadium billiard. Our results are twofold. First, we 
develop a method to solve in a simple way the quantum mechanical evolution 
of planar billiards with moving walls. The formalism is based on the scaling 
method [I] which enables the resolution of the problem in terms of quantities 
defined over the boundary of the billiard. 

The second result is related to the quantum aspects of dissipation in systems 
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with complex spectra. We conclude that in a slowly varying evolution the e- 
nergy is transferred from the boundary to the particle through Landau— Zener 
transitions. 
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INTRODUCTION 



The way in which energy is transferred from the time dependent mean field to the indi- 
vidual nucleons is an important ingredient, for instance, in fission processes |2|] and in large 
amplitude collective motion at low energies. Because of the Pauli principle it is expected 
that one body effects, i.e., loss of energy due to collision of independent individual nucleons 
with the mean field, should dominate the dissipation mechanism. 

Several descriptions of these processes involving different approximations are available in 
the literature |3|-|5[]. These theories are perturbative in character and linear in the collective 
motion. Therefore they are not suited to address the issues related to nonlinear dynamics 
and the onset of chaos. On the other hand, the integrable or chaotic nature of the motion 
is of crucial importance to the dissipation mechanism |7j],z.e., the transition from order to 
chaos provides the possibility for a variety of nuclear responses (from elastic to elastoplastic 
to dissipative). Therefore detailed studies of simplified models resembling a nucleus may be 
of interest. 

Planar billiards are perhaps the best systems to model the processes described above in 
which the nucleus can be imagined as a time-dependent container filled with a gas of non in- 
teracting point particles 0|| . Billiard systems have been thoroughly studied in the context 
of classical and quantum chaos ||. In particular it has been shown that the quantum spectra 
of generic planar billiards have GOE (Gaussian Orthogonal Ensemble) characteristics that 



are observed in the excited spectra of nuclei |10 



In a seminal paper, Hill and Wheeler |TTJ suggested the Landau— Zener (L-Z) transitions 
as a mechanism for nuclear dissipation. The mechanism is based on the excitation of the 
individual nucleons via transitions at avoided level crossings near the Fermi surface. These 
excitations produce the damping of collective motion describing deformation of the nucleus. 
In an adiabatic evolution of the collective coordinates, the nucleus changes its shape rela- 
tively slowly, while the nucleonic levels move up and down in energy. Small deformations 
in the nuclear shape produce occasionally that two nucleonic levels almost cross each other 
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and experience an avoided crossing. During the whole process many avoided crossings occur 
with more or less random transitions between nearest neigbours, in such a way that the 
system may end up in an arbitrary energy state. Within this picture one can imagine a 
stochastic dynamics in which by simply reversing the temporal evolution one does not re- 
cover the initial state. Therefore, the internal degrees of freedom are excited and the motion 
of the collective coordinates is thus damped. 

More recently Wilkinson, making good use of the properties of complex spectra, introduced 
a statistical treatment of dissipation in finite-sized quantum systems in terms of L-Z tran- 



sitions in the context of random matrix theory [12|. However, a L-Z mechanism as the 



generator of dissipation in systems with complex spectra, has recently been seriously ques- 



tioned [ 13 1 . In their works Bulgac and collaborators suggest that the diffusive process in 



energy is dominated by memory effects and that the picture for dissipation through L-Z 



transitions is likely to be incorrect (see section [III] for a detailed analysis of these argu- 
ments) . 

Obviously the best way to elucidate this question is to solve the quantum mechanical evolu- 
tion of a generic system, though this is difficult even for planar billiards with moving walls. 

The goal of the present article is to present a new formulation to solve the quantum 
mechanical evolution of planar billiards with moving boundaries. Using this formulation 
we study the time evolution of a specific billiard, the Bunimovich Stadium with externally 
driven walls. We restrict the analysis to slowly varying (adiabatic) evolutions ||14|| . The 



notion of slow motion will be quantified in section [TV|. Our aim is to understand whether 
a L-Z mechanism generates the damping of the slow degree of freedom in a system with 
complex spectra. 

The work outline is as follows. In the next section we introduce a one-dimensional new 
formulation in order to solve the quantum evolution of planar billiards with moving walls . 
Section [D] is devoted to the numerical results obtained for the Bunimovich Stadium billiard 



with GOE spectrum characteristics. Using these results, we evaluate relevant properties for 
the dissipation process. In section [TTT] we discuss in detail whether a L-Z transition mecha- 
nism describes the dissipation of the slow degree of freedom or equivalently the difussion of 
the fast ones. The last section is devoted to final remarks and conclusions. 
Before proceeding we want to stress that planar billiards systems externally driven can also 
be used to model other problems often encountered for example in mesoscopic systems, 



atomic clusters and of course deformable cavities ||15|| . Therefore, the results presented in 



this work could also contribute to domains other than nuclear physics. 

I. THE METHOD 

In a recent work Vergini and Saraceno developed a method to calculate directly all 
eigenvalues and eigenfunctions in a narrow energy range of quite general time independent 
2 — d billiards, by solving a generalized eigenvalue problem in terms of quantities defined 
over the boundary. The method is based on the use of scaling that enables to write the 
boundary norm explicitly as a function of the energy 

The aim of this section is to extend the method of scaling to solve the Schrodinger equation 
for 2 — d billiards with time dependent boundary conditions. 

Let C(t) be a closed curve defining at time t a two dimensional domain T>(t). We restrict to 
star shaped domains, this means that r n = r.n > Vr G C(t)\ n is the outgoing normal to 
C(t). Consider a particle of mass m inside the billiard, then the Schrodinger equation reads, 

?f = l -A*. (1.1) 
dt 2m K ' 

^ satisfies the time dependent boundary condition \I/(C,£) = where ( is a point on C(t), 
and we consider functions normalized to one on the domain. A standard procedure is to 
expand the solution in terms of the adiabatic basis, 

*(r,*) = £cv(t)P M (t)<^M). (1.2) 



5 



Pfj,(t) = exp(— i ^uj^dt') is the contribution of the dynamical phase with to^ = %k 2 (t)/ 2m. 
The adiabatic basis {0^} constitutes a complete set of real eigenfunctions of the billiard 
at each time; that is, cf)^ satisfies the Helmholtz equation A0 /t (r,t) = — /c^(t)0 A1 (r, t) with 
Dirichlet boundary condition M (£,t) = 0, and it is a continuous function of time. 
We generate from the original domain defined by C(t) a family of systems that depends on 
a parameter a. These systems evolve with the curves C a (t) that are obtained from C(t) 
through a scaling transformation on the plane r — > ar (if ( is a point on C(t), then (/a is 
the corresponding point on C a {t)). 

To each function M (r, t) we associate the scaling function M (a, r, t) = M (ar, t). This family 
of functions depending on the scaling parameter a verifies the Helmholtz equation with wave 
number ak^ and satisfies the Dirichlet condition over the scaled boundary. Moreover, we 
require that the mass of the particle in the scaled systems changes as a 2 m in such a way 
that uj^ is independent of a. The last statement implies that the time evolution is the same 
for all the scaled systems. 

Our approach to solve the Schrodinger equation is to expand the wave function in terms of 
the adiabatic basis represented by the scaling functions. After replacing the expansion into 
the equation ( |1.1| ) we obtain, 

5>„(t) P v (t) <i> v {a,T,t) = - J2 a »(t) P»(t) ^(a,r,f) . (1.3) 
Differentiating this equation with respect to a results in, 

P v {t) ^(<*,r,t) = -D^i) P„(t) f^(«,M) • (1-4) 

The remainder of the calculus consists on commuting the order of the partial derivation 
in the rhs of (|1.4j) , multiply the equation by dcp^/da (a,r,t) and specialize the resulting 
equation in a = 1. Finally we integrate over the boundary of the billiard C(t). After this 
straightforward calculation, the final equation reads, 

t+\ r> u\ I dfadfodl ^ r <% d d(j> v dl 

V auit) m U = ~V aAt) W) h) ~ n ' (L5) 



where dl is the length element on the boundary. For the sake of simplicity we have omitted 
the argument (a = l,r, t) in the last equation. In |]J it was proved that the integral in the 
lhs of the last equation verifies a quasiorthogonality relation, this means 



bp d(j) u dl (k^ k v ) 



2kl Jcit) da 8a r n 6ftv + (k, + ' (1 ' 6) 

Employing this important relation in ( |1.5|) , we derive the standard system of differential 
equations in the adiabatic basis 

^(t) = -Y,<t){Pv{t)/PM c^t) (1.7) 

V 

with the coefficients C^ u computed approximately in terms of quantities defined over the 
boundary, 

W^wjw-toditt^- (L8) 



The exact expression for the coefficients follows from (|1.3|), 



C» v {t)= <f>,{v,t) ^(r,t) da . (1.9) 
Jv{t) ot 

In the last equation each C^t) involves an integration on the domain T>(t). It is also very 
easy to prove that they are antisymmetric. 



To compute each via the equation (|1.9|) the domain of integration has to be partitioned 
at least in iV pa k 2 cells, with k equal to the maximum among the wave numbers of the func- 
tions in the region of integration, k = max{k^, k u }. If n is the dimension of the adiabatic 
basis restricted to the energy region where the evolution will take place, one needs to know 
n 2 coefficients. As a consequence, the dimension of the problem of finding the coefficients 

from equation ( |1.9| ) turns to be of the order of n 2 x iV at each time. 
As a way to check the goodness of the new formulation we have computed for the specific 
billiard studied in section [TI|, the coefficients C^it = 0) for a fixed fi, with v = fi + j 
(j = 0, ±1,±2, ...) calculated exactly (equation (|1.9|)) and using the equation ( |1.8| ) (see 
Fig. |l|). The correspondance is extremely good over a great number of levels. The departure 



between the two plots begins for | j |~ 10, but in this region the values of the coefficients 
are very small. 

With the present formulation the CPU time necessary to compute the coefficients is con- 
siderably reduced in comparison with the time needed in the standard approach (equation 
( |1.9[ )). From the preceding remarks, and in order to study the interaction between neigh- 
bouring levels in the spectrum, we will calculate the coefficients employing the relation 

(ED- 

II. NUMERICAL RESULTS 

Using the method presented in section |, we will analize the dynamics of a particle of 
mass m inside a Bunimovich stadium billiard with moving boundaries. 
A point particle inside the static stadium billiard is a very well known example of a fully 
classical chaotic system |20[ . The particle moves freely on the two dimensional domain and is 
perfectly reflected from its boundary. The boundary is formed by two semi-circles of radius 
r connected by two straight lines of length 2a. Fig. |^ shows a desymmetrized version of the 
system with area 1 + 7r/4. 

To study the dynamics, the parameter i = a/r is changed with a finite velocity I in such a 
way that the total area of the billiard remains unchanged. We have fixed the area to avoid 
a drift in the energy spectrum; this situation is characteristic of nuclear processes where the 
nucleonic density is approximately constant. The drift term represents a reversible change 
in the energy of the system and can be neglected in the analysis of an irreversible dissipation 



process [[H]]. Therefore, the dynamics of the boundary is introduced through the function 

e(t). 



Fig. |3] shows the spectrum of k = y 2mE / h as a function of £, 1 < I < 1.14. We have 
selected the wave numbers k^i) between 48.8 and 50 because in this region a large number 
of energy levels exists in a narrow portion of the spectrum. Although the properties that we 
are going to evaluate are characteristic of this region of the spectrum {k ~ 50), as we will 



S 



show below, a proper scaling can be done in order to evaluate them in other energy regions. 
The spectrum exhibits the typical behavior of avoided levels crossings that characterizes the 
energy levels as a function of a parameter for general systems without constants of motion 



22|| . Also, we recognize that some avoided crossings are situated on two parallel lines labeled 



L3 and L4 (see Fig. |3p. These lines are associated to bouncing ball states with three and four 



low excitations respectively. These states are highly localized in the momenta space |f23 



therefore their interaction with neighbouring states is smaller than the interaction between 
delocalized generic states. 



Let us analize the coefficients C^ v which determine the quantum mechanical evolution of 
the system (see section |). They may be expressed in terms of the deformation parameter 
£, as they satisfy : C^ u {t) = iC^i). 

Fig. |] shows the functions | C w+ i(£) | for several pairs of nearest neighbouring levels. A 
well defined structure of peaks is observed. The peaks appear each time two neighbouring 
energy levels experience an avoided crossing (it is very easy to follow in the Fig. |3] a pair 
of energy levels as a function of the parameter I in order to confirm this assertion). The 
height of the peaks diminishes when the energy gap between levels at the avoided crossing 
increases. For this reason the peaks corresponding to interaction with bouncing ball states 
are one order of magnitude greater than the generic ones. We label with a, b, c, small 
peaks that correspond to not well defined avoided crossings or to situations where it is still 
difficult to decide whether an avoided crossing exists by simple inspection of the spectrum 
(see also Fig. [3D. 

For second neighbouring levels, we also find some well defined peaks, they appear essentially 
when three levels come close to each other (this situation is discussed in the next section). 
The heights of the C wt+2 (^) peaks are one order of magnitude smaller than those of the 

For coefficients with |/i — v\ > 2, we do not observe any simple structure, however the am- 
plitude of these coefficients is indeed very small, lower than five in the scale of Fig. [I], in 
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comparison with the amplitude observed for nearest neighbouring coefficients. 
From the present analysis it is clear that the information contained in the coefficients en- 
ables a complete definition of the avoided crossings and that this information is not always 
available in the spectrum. 

The peaks between first neighbouring levels are very well fitted by lorentzian functions, 
as it is expected for a L-Z transition (the first part of the next section is devoted to explain 
the expected lorentzian behavior of the coefficient), for almost all peaks with the exception 
of some small ones. 

Fig. |5] summarizes the preceding remarks. It shows on the top the function | C13 12 (-0 | 
—Lz(£) where Lz(£) is a sum of lorenztian functions centered on the well defined peaks 
of the coefficient Ciz^t)- Each lorentzian function is defined in the next section by the 
equation ( |3.1| ). The widths and the position of the centers are £i nt and £0 respectively. 
The remainder of the Fig. [| shows the functions | Ci 3 vi-j{£) I (j = 2, 3, 4) as a function of 
£. The figure reveals the lorentzian behavior of the first neighbouring levels coefficients, and 
the lack of a defined structure in the coefficients C^(£) for |/i — u\ > 2. 

The previous numerical study would be still more appealing if we knew how the spec- 
trum scales to other energy regions. The Weyl's law |24| tells us that the density of states 
associated to the vertical axis in Fig. |3|, scales as k. The problem appears with the horizontal 
axis because the scaling of the density of consecutive avoided crossings p a c is unknown. 
Working in different energy regions and after an exhaustive numerical analysis, we have 
obtained that p a , c scales as k d with d = 1.92 ± 0.1. Fig. |6] shows p a . c as a function of k for 
gap sizes less than one quarter of the mean level spacing. In this calculation we have not 
considered avoided crossings with bouncing ball states because their relative contribution to 
the density of states decreases as k^ 1 ^ 2 |f25f . 

Another important fact to stress is that each peak is very well defined; its width (given by 
£i n t) is much smaller than the mean distance between consecutive peaks Pa.c" 1 - Only in few 
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cases a small overlap between consecutive avoided crossings is observed (see for example 
the peaks labeled by in Fig. For generic avoided crossings we have obtained that 
Pint pa.c ~ 0.2; for avoided crossings with bouncing ball states this product is even smaller 
as expected. 



III. LANDAU- ZENER BEHAVIOR 

As we have mentioned in the introduction, our aim is to understand whether L-Z tran- 
sitions govern the mechanism of one body dissipation in systems with complex spectra. 
In the previous section we have analyzed the coefficients C^ v that describe the quantum 
evolution of a particular system with complex spectrum. The analysis revealed that the 
coefficients have a simple structure of well defined lorentzian peaks as the dominant contri- 
bution, plus a very small component without any defined structure (see Fig. [5]). These peaks 
are concentrated in the first neighbouring levels coefficients, which, as we will show below, 
is a characteristic of L-Z transitions. In few cases well defined peaks appear in coefficients 
between second neighboring levels, but these peaks are one order of magnitude smaller than 
the previous ones. We will discuss this situation below using an idealized three level system. 



We begin this section with a brief review on the theory of L-Z transitions. Consider 
a two level system that depends on a parameter £, in such a way that for £ = £q the 
energy levels experience an avoided crossing. Let 4>±{£) be the adiabatic eigenstates and 



E + — E- = \J^ 2 {£ — £q) 2 + 4e 2 the energy gap between the associated eigenvalues, with 7 
and e constants. The adiabatic theorem |18|] tells us that if the system is initially in the state 
0_ and £ changes infinitely slowly from £ < £q to £ > £q the system will remain in the state 
0_. However, if £ changes with a finite velocity the final state will be a linear combination 
of the basis states. Zener derived the probability of an adiabatic transition employing the 
diabatic basis |T7| for a constant velocity of the parameter £. If at time t = —00 the system 
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were in the state 0_ the transition probability at time t = oo is P z = exp{—2T[e 2 /^j£K) JT 
Using the adiabatic basis and ( |1.9|) is straightforward to derive 



where li nt = 2e/ r y is the width of the lorentzian function; that is, the characteristic time for 
a L-Z process. 

It has been suggested in recent literature that it is very difficult to characterize the inter- 
action between neighbouring levels in spectra like the present one in terms of L-Z transitions 
13| . The arguments could be summarized as follows: i) It is not always possible from the 



spectrum to localize the position of the avoided crossings and to determine the parameters 
that define the L-Z transitions. This assertion is partially true; it is not in the spectrum 
where all the information is contained. We have solved this problem employing the adiabatic 
basis, in which the position of the avoided crossing and the interaction length £ int are well 
defined in terms of the coefficients C^ u . 

ii) The L-Z transition probability is exponentially small when the length of the transition 
process goes to infinity, but it could be strongly affected for lengths of the order of £ int [P5 



This problem could emerge if the mean distance between avoided crossings p a , c is of the order 
or less than £ int . In the preceding section we have obtained £ int p a . c ~ 0.2 for generic avoided 
crossings (between delocalized eigenfunctions) and this value is highly reduced for localized 
eigenfunctions. In a physical system the eigenfunctions present some degree of localization 
because the associated classical phase space is not fully chaotic. Therefore, we do not expect 
correlations between consecutive avoided crossings. In terms of the coefficients C^+i, the 
last assertion means that each individual peak is very well defined. 

iii) One often encounters a situation where three levels come close to each other and by 
simple inspection one can think in a three level crossing (see for example points A, B, C 
and D in the Fig. |3|). In order to understand this process we will analyze a three level system 
which mimics such a circumstance. 
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Consider a one parameter dependent Hamiltonian defined in the diabatic basis by the fol- 
lowing matrix, 

/ ' '"int 1 ^ 



1 1 

^ 1 £/£ int J 

with e the perturbation and £i n t the characteristic transition length. This Hamiltonian can be 
diagonalized analytically for each time. The upper and lower eigenenergies are represented 



by hyperbolas E± = ±€^(£/£ int ) 2 + 2 , as in the L-Z process, and the middle energy is 
E = for all times (see Fig. [7]). Obviously, for a diabatic evolution, if the system were 
in the upper stcite dbt — ^l^int 

» 1, there is a high probability ~ 1 that the system decays 
to the state £L at £/£i n t » 1. In other words, the presence of £L affects enormously 
the transition probability between E + and Eq. However, we are interested in an adiabatic 
evolution where the transition probability to E_ turns to be small. In such a situation we 
want to determine whether the L-Z parameters for E + and Eq, that is (E + — E )/2 at £ = 
(we denote it A) and £i nt , adequately describe the transition probability between these two 
states. This point is not obvious at all. For example, the distance (E + — Eq) is largely 
affected by the presence of E- and A/e, contrary to the case of a two level system. 
We have computed numerically the dynamical evolution of this model obtaining the following 
result for the transition probability between E + and Eq: 

Pe + ^e ~ \ exp[-vr A {£ int / IK) 0.96] (3.2) 

for Pe + ^e ~ 0.2. That is, although the parameters need to be renormalized, the factor is 
very close to 1. As a conclusion, this three levels system may be thought as two independent 
avoided crossings with L-Z interactions. 



IV. FINAL REMARKS 



The results of this paper attempt to extend the present understanding of one body dissi- 
pation processes. To analyze the way in which energy is transferred form the time dependent 
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mean field to the individual nucleons we have modeled the mean field by a slowly time de- 
pendent container. The same approach has been employed by many other authors in related 
nuclear models ]?||| . The container is represented by a planar billiard with externally driven 
moving walls. To solve the quantum mechanical evolution of these simplified systems, we 
have derived a one-dimensional formulation. This approach gives the possibility to study 
the evolution of highly excited states and reduces the CPU time involved in the calculations. 
We have devoted part of the work to answer a fundamental question; whether a 
Landau— Zener excitation mechanism governs the irreversible transport of energy from the 
driven wall to the particles in an adiabatic evolution. We have analyzed a parameter depen- 
dent billiard system with GOE character spectrum, concluding that in an adiabatic evolution 
of the external parameter, the dissipation is dominated by L-Z transitions at the avoided 
crossings. The adiabatic limit is attained in the limit of an infinitely slow evolution. On 
the other hand, adiabatic evolution refers to slowly varying evolutions |14|. Of course, the 
notion of slow motion needs to be clarified. For example, we have excluded in our analysis 
the structure showed by the function Ci 3 n(£) in Fig. |5] because its height is very small; 
although the area under it is comparable to the area under any peak observed in Fig. |4j. 
However, because Ci 3 n (£) is multiplied in the differential equation fll.Tj ) by an oscillatory 
function with period T « hpE {pE is the density of energy levels), its effective contribution 
is canceled if the time required by the collective motion to sweep the structure t co u ~ p~^}J£ 
is larger than T. 

The above adiabaticity condition is satisfied by systems where quantum effects are very 
important, such as nuclei. However, as the wave number increases, the collective veloc- 
ity needs to be reduced drastically. Taking the semiclassical limit H — > 0, k — > oo, with 
(kh) 2 /2m = E = const , it results that T = O(k ) and t coll = 0(k~ d /l) . Therefore, for 
any finite value of £, the evolution is always diabatic in the semiclassical limit. In another 
terms, a semiclassical theory of dissipation requieres a scaling of £. To our knowledge, this 
important point has not been taken into account in previous works W3 . 
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The description of the damping process in terms of L-Z transitions has been already done 
by Wilkinson in the context of pure random matrix theory [P^,2Tf| but to our knowledge the 
present work is the first study carried out for a more realistic system. If the L-Z behavior 
holds it is more or less straigthforward to write the corresponding diffusion equation to 
quantify the dissipation mechanism (for further details see 
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FIG. 1. Coefficients (with = 49.456279) as a function of k v computed exactly (equation 
(4)) and employing the boundary definition (equation (8)). The system used is introduced in 
section |Il| and the calculation corresponds to t = 0. 
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FIG. 2. Desymmetrized Bunimovich stadium billiard. The area of the billiard is fixed to the 
value 1 + 7r/4. Then, the boundary only depends on one parameter (£ = a/r). 
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FIG. 3. Spectrum of the Bunimovich stadium billiard as a function of £, 1 < I < 1.14. The 
wave numbers k^(t) run between 48.8 and 50. See text for more details. 
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FIG. 4. | C M ^ + i | as a function of i for several energy first neighbouring levels. The labels a, b, 
c, show small peaks that correspond to avoided crossings whose parameters can not be obtained 
directly from the spectrum. The peaks that correspond to interaction with bouncing ball states are 
out of scale and their maximum values are shown. The labels A, B, C and D show well defined 
peaks that in the spectrum appears as three levels avoided crossing. The label O shows few cases 
where there is some overlap between consecutive avoided crossings. 
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FIG. 5. The top shows | C13 12 {t) \ —Lz(£) as a function I. Lz(£) is a sum of lorenztian functions 
centered on the peaks of the coefficient | C13 |. Their widths and the position of the centers 
are 4nt an< i ^0 respectively. The remainder part of the figure shows the coefficients | C13 i3-j(£) \ 
for j = 2, 3, 4 as a function of I. 
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FIG. 6. Log-Log plot of density of avoided crossings p a . c as a function of k for gap sizes less 
than one quarter of the mean level spacing. p a . c scales as k d with d = 1.92 ± 0.1. 
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FIG. 7. Energy levels E+, E- and E$ as a function of ljl%nt for the 3— levels hamiltonian 
analized in section p| . The dotted lines indicates the asymptotes to E— and E+. The distance 
{E + — Eo)(0) = \/2e is also drawn. 
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